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Abstract. We prove certain duality properties and present recurrence relations 
for a four-parameter family of self-dual Koornwinder-Macdonald polynomials. The 
recurrence relations are used to verify Macdonald's normalization conjectures for 
these polynomials. 



I. Introduction 

In a to date unpublished but well-known manuscript, Macdonald introduced certain 
families of multivariable orthogonal polynomials associated with (admissible pairs of 
integral) root systems and conjectured the values of the normalization constants 



turning these polynomials into an orthonormal system [Ml]. Recently, Cherednik 
succeeded in verifying Macdonald's normalization conjectures in the case of reduced 
root systems (and admissible pairs of the form (R,R V )) using a technique involving 



so-called shift operators [|C 1|| . Previously, this same technique had enabled Opdam to 



prove the normalization conjectures for a degenerate case (q — > 1) of the Macdonald 
polynomials known as the Heckman-Opdam-Jacobi polynomials |0|, |H|]. 

Meanwhile, a generalization of Macdonald's construction for the nonreduced root 
system BC n — resulting in a multivariable version of the famous Askey- Wilson poly- 
nomials [ |AW| ] — was presented by Koornwinder [ K2[| . It turns out that all Macdonald 



polynomials associated with classical (i.e., non-exceptional) root systems may be seen 
as special cases of these multivariable Askey- Wilson polynomials |D1[ , Sec. 5] (type 
A by picking the highest-degree homogeneous parts of the polynomials and types B, 
C, D, and BC, by specialization of the parameters). 

In the present paper, we will prove certain duality properties and recurrence rela- 
tions for (a four-parameter subfamily of) the Koornwinder-Macdonald multivariable 
Askey- Wilson polynomials, which enable one to verify the corresponding Macdonald 
conjectures for the (ortho)normalization constants also in this (more general) situa- 
tion. Our approach does not involve shift operators but rather exploits the fact that 
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the polynomials are joint eigenfunctions of a family of commuting difference opera- 
tors that was introduced by the author in Ref. |pi|| (see also Ref. |p3| ) . By duality, 
these difference operators give rise to a system of recurrence relations from which, in 
turn, the normalization constants follow. 

The same method employed here was used already several years ago by Koorn- 
winder when verifying similar duality properties and normalization constants for the 
Macdonald polynomials related to the root system A n \ K1 , |M3[| . (In this special 
case, though, the validity of the normalization conjectures had also been checked by 
Macdonald himself.) The A n -type Macdonald polynomials constitute a multivariable 
generalization of the g-ultraspherical polynomials | ]AW| | (to which they reduce for 
n = 1). The present paper may thus be regarded as an extension of Koornwinder's 
methods in Ref. [ Kl|| (see also Ref. [M3, Ch. 6]) to the multivariable Askey- Wilson 
level, or, if one prefers, as an extension from type A root systems to type BC root 
systems. 



2. Koornwinder-Macdonald polynomials 

The Koornwinder-Macdonald multivariable Askey- Wilson polynomials are charac- 
terized by a weight function of the form 

A(x) = I [ d v (exj + e'xf) j [ d w (exj), (2.1) 



where 



l<j<j'<n ±<j<n 
e,e'=±l e=±l 



dJz) = - — , q = e p , 

(qBe- a ",q) O0 ' 

d (z) = {e~ az , -e~ az , q x l 2 e- a \ -q^e-™; g)^ 

w (q9°e~ az , —q 91 e~ az , q(9t +1 / 2 )e~ az , — g(9 3 + 1 / 2 )e~ a2; ; q) 



oo 



and the g-shifted factorials are defined, as usual, by (a;g)oo = n^o(-'- — a( l l ) anc ^ 
(ax,... ,ak,q)oo = {a-i', q)oo • • • (ofe! o)oo- To ensure the convergence of the infinite 
products contained in A (|2.1|) it will be assumed that a, (3 > (so 1 < q < 1); in 
addition, we will also assume g, g r > 0, r = 0, 1, 2, 3. 

Let {m\(x)}\(zA denote the basis consisting of even and permutation symmetric 
exponential monomials (or monomial symmetric functions) 

m x (x) = e ° E ' ,=1 X ' iXj > A G A = {A G Z n | Ai > A 2 > ■ • • > A n > }, 

y&wx (2.2) 

with W being the group generated by permutations and sign flips of Xj, j — 1, . . . ,n 
(W — S n x (Z2) n ). The monomial basis can be partially ordered by defining for 
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A, A' e A (O 



A'<A iff Yl X 'i ^ Yl X i for ™ = l,...,n (2.3) 

l<j<m 

(and A' < A iff A' < A and A' 7^ A). The Koornwinder-Macdonald polynomials p\{x), 
A G A can now be introduced as the (unique) trigonometric polynomials satisfying 

i. p\(x) = m\(x) + ^2 c a,a' mx'(x), c a ,a' e C; 

A'eA,A'<A 

ii. (px,m X ')A = if A' < A, 

where (•, -)a denotes the inner product determined by 

(m x , 771a')a = ( 77- ) / ■••/ m\(ix)m\i(ix)A(ix)dxf-dx n 

\Ali/ J_ n / a J-ir/a ^2.4) 

(and extended by bilinearity) . In other words, the polynomial p\(ix) consists of the 
monomial m\(ix) minus its orthogonal projection in L 2 (] — -] n , A(ix)dxi ■ ■ ■ dx n ) 
onto the finite-dimensional subspace span{mA'(ix)}A'eA, a'<a- 

It is of course possible to extend this construction of multivariable polynomials 
determined by Conditions i. and ii. to a more general class of weight functions than 
the one considered here. In general, however, the resulting polynomials will not be 
orthogonal (except for n = 1) because the ordering in Eq. fl2.3|) is not a total ordering 
(unless n — 1). (A priori the construction only guarantees that p\(x) and p\>(x) be 
orthogonal if A and A' are comparable with respect to the ordering in Eq. (|2.3|).) Still, 



it turns out |[K2|| that for the weight function A (|2.1|) the corresponding polynomials 
indeed do constitute an orthogonal system for arbitrary n: 
Orthogonality 

(Pa,Pa')a = if A^A'. (2.5) 

This feature should be looked upon as a very restrictive property of the weight func- 
tion a (grp. 

Remark: The relation between our parameters and the parameters employed by 
Koornwinder reads (cf. ||K2| , Eqs. (5.1), (5.2)]) 

t = q 9 , a = q 90 , b = -q 91 , c = q i92+1/2 \ d = -q {93+1/2) . (2.6) 

Furthermore, Koornwinder fixes the period the trigonometric functions to be 2ir(i), 
i.e., he puts a = 1. 



3. Difference equations 



Another special property of the polynomials associated with A (|2.1| ) is that they 
satisfy a second order difference equation |[K2|| . (This difference equation is in fact 
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instrumental in the orthogonality proof.) It can be written as 



with 



and 



l<j<n k^j 

w {~ x i) \\ v (~ x j + x k) v(-Xj - x k ) \p x (x - (3ej) - Px(x)\\ 

= E(p + X) p x (x), 



v(z) 



sh§0% + z) 



w(z) 



sh(fz) 

shf (f3g + z) chf (/%! + z) shf (/9</ 2 + § + z) chf (/fy 3 + f + z) 



sh(fz) ch(fz) shf(f + z) chf(f + z 



£(y)=2 (ch(a/%)-ch(a/%)), 
i<i<« 

P= X] Pj e 3> Pj = ( n -j)9+(9o + 9i + 92 + 93)/^- 

l<j<n 



(The vector ej denotes the j-th unit element of standard basis in K n .) For n = 1 
this difference equation reduces to the well-known difference equation for the Askey- 
Wilson polynomials |AW|. 



In Ref. |pi|| it was shown that for arbitrary number of variables n the above 
difference equation can be extended to a system of n independent difference equations 
of order 2r, r = 1, . . . , n, respectively. This system is explicitly given by: 

Difference equations 



■ ,"},0<|J|<r 
=±l,j£j 



Ujc jr _\j\(x) V £j>J c(x)p x (x + f3e £j ) = E r (p + \)p x (x), r 



(3.2) 
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with 



V £j , K (x) = Y[w(ejXj) Y[ v(ejXj + eyx^vfaxj + e r x r + (5) 

i<i' 

Y[ v{ejXj + x k )v(ejXj - x k ), 



j<j' 



Uk, p ( x ) = X] Ii^C^O 1 [ v(£iXi + e V X V )v{-EiXi - EyXy - (3) 



LcK,\L\=p leL IJ'eL 
ej=±l, leL 1<V 



I f v(eixi + x k )v{eixi - 



Xk) 



leL 

k&K\L 



and 



E r (y) = 2 r (-l) rHJ| (n ch M%') ch(a/3p, 1 )...ch(a^ lr _ |J| 



Jc{l,— ,n} jeJ r<h<—<l r -\j\<n 

0<\J\<r 

In the above formulas \J\ represents the number of elements of Jc {1, . . . , n} and 

e £j = J2 £ j e 3 (£j = ±l)- ( 3 - 3 ) 

Furthermore, we used the conventions that empty products are equal to one, Uk, p = 1 
if p — 0, and the second sum in E r (y) is equal to one if \J\ = r. For r = 1, the 
difference equation in Eq. (|3.2| ) reduces to that of Eq. (|3.1| ). 

Remarks: i. Equation ( |3.2| ) may be interpreted as a system of eigenvalue equations 



(D rPx )(x) = E r (p + A) Px (x), (3.4) 
for a family of n independent commuting difference operators of the form 

D r (x)= ^2 Ujc^\j\(x)V £ j t jc(x)T £ j tPl r = l,...,n, (3.5) 

JC{1,... ,n},0<| J|<r 
£,=±1, jeJ 

with T sJi/3 = // ,.; and 

(T±j t /3f)(xi, . . . , x n ) = /(xi, . . . , ± (3, Xj + i, . . . , x n ). 

The commutativity of Di(x), . . . ,D n (x) follows [[L>1|| from the fact that the differ- 



ence operators are simultaneously diagonalized by the basis {p\(x)}\ e A- Two other 
properties of the operators that follow immediately from their diagonalization by the 
Koornwinder-Macdonald polynomials are the triangularity with respect to the par- 
tially ordered monomial basis {m\(x)}\(zA, and the symmetry with respect to the 
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inner product (•, -)a (H 
has (Triangularity) 



(the eigenvalues E r .(p + A) are real). More precisely, one 



(D r m x )(x) 
and (Symmetry) 



E r (p + A) m x (x) + ^ [D r ]x,x'm y (x), [D r ] x ,y E C, 



x'eA,\'<\ 



(D r m x ,m X i)A = (m x ,D r m y ) A . 
ii. For r = 1 the operator D r (x) (|3.5| ) reduces to (cf. the l.h.s. of Eq. (|3.1| )) 

D ( x ) = Yl ( w ( x i) II V ( X J + v ^ x i ~ Xfc ) ( T i>P ~~ X ) + 



(3.6) 
(3.7) 

(3.8) 



r 



l<j<n k^j 

w(-Xj)Y[v(-Xj + x k ) v(-Xj - x k ) (T_ jt p 

With the aid of this operator the following useful representation for the Koornwinder- 
Macdonald polynomials can be given (cf. Ref. |M1|| ) 

D(x)-E(p + X / ) 



P\{x) 



n 



E(p + A) - E(p + A') 



m x {x) 



(3.9) 



a'ga,a'<a 

(Using the Triangularity (|3~6| ) and Symmetry ( |3.7|) of D(x) (3.8), it is not difficult 
to verify that the r.h.s. of Eq. ([3.9)) satisfies the Conditions i. and ii. in Section ||. 
Notice also that the denominators in the r.h.s. of Eq. (|3.9|) do not vanish because 



E(p + A') < E(p + A) if A' < A (cf. Ref. [0, Lemma 5.1]).) It follows from this 
representation for p x (x) that the Koornwinder-Macdonald polynomials are rational 
in exp(aj3g) and exp(a/3g r ), r = 0, 1, 2, 3. Hence, they may be extended uniquely to 
nonnegative (or even complex) values of the parameters g, g r . In view of the analytic 
dependence on the parameters it is clear that the resulting polynomials then satisfy 
the Difference equations (|3.2|) for all the values of g, g r , r = 0, 1, 2, 3. 



4. Duality and recurrence relations 

Note: In this section, we will drop the condition that the parameters g, g r , r = 
0, 1, 2, 3 be nonnegative (cf. Remark ii. of the previous section). 

In order to describe the duality relations it is convenient to introduce certain dual 
polynomials p x (x), A G A. These dual polynomials are again Koornwinder-Macdonald 
polynomials but with a slightly different parametrization. Specifically, the parameters 
of p* x (x) are related to those of p x (x) by 



9* 



a, 
9, 



[9o\ 




/ 1 


1 


1 


1 \ 


(go\ 




1 


1 


1 


-1 


-1 


9i 


9*2 


~ 2 


1 


-1 


1 


-1 


92 


w 




V 1 


-1 


-1 


l ) 


w 



(4.1) 
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Notice that the reparametrization in Eq. ( |4.1| ) is involutive, i.e., p*\{x) = p\(x). 

Furthermore, instead of working with monic polynomials we go over to a different 
normalization by introducing 

(p* = YTj=i Pj e j with p* = {n — j)g* + (flft + ' — 1~ ^3) /2) , where we have also rescaled 
the arguments of the trigonometric polynomials such that p\(x) and p\{x) have the 
same (imaginary) period in Xj, j = 1, . . . , n (viz. 2iri/(a(3)). Of course, the renor- 
malization in Eq. ( |4.2| ) only makes sense provided P\{f3p*) and p\{dp) do not vanish. 
This is guaranteed, at least for generic parameters, by the following lemma. 

Lemma 4.1. For generic parameters one has 

p x ((3p*),P* x (ap)^0. 

Proof. For g,g , . . . ,g^ = 0, the polynomial p\(x) reduces to the monomial symmetric 
function m\(x) and p* = 0. Thus, it is clear that for this special choice of the 
parameters p\{(3p*) 7^ 0. But then the same follows for g,g , . . . , g 3 in an open dense 
subset of R (or C) because of the analytic dependence of p\(x) on the parameters (cf. 
Remark ii. in Section [3]). The analogous statement for p* x (ap) follows by duality. □ 

In Section [5] the value of p\{j3p*) will be computed explicitly. We will then see that 
px(Pp*) is positive for all nonnegative values of parameters g, g r , r = 0, 1, 2, 3. 

The matrix in Eq. ( |4.1| ) relating (g$, . . . ,g^Y and (go, ■ ■ ■ , g^Y has eigenvalues +1 
(with multiplicity three) and —1 (with multiplicity one). The invariant subspace 
corresponding to the eigenvalue +1 consists of the hyperplane g — g\ — g 2 — g-s = 0. 
For parameters in this hyperplane one has g* = g r (r = 0, 1, 2, 3) and p\(x) = p\(x). 
In other words, for these parameters the polynomials p\(x) are self-dual. In the rest 
of the paper we will always assume that the Self-duality condition 

9o-9i- 92 - 93 = (4.3) 

is satisfied (unless explicitly stated otherwise). 

After these preparations we are now ready to formulate the duality theorem, which 
relates the value of fi\(x) in the point p* + p to value of p^x) in the point p + A 
(A, p G A). 

Theorem 4.2 (duality relations). Let A, p G A Q2.2p . Then the renormalized Koorn- 
winder-Macdonald polynomials p\(x) and p*^(x) ( |4.2| ) satisfy the relation 

p x (p* + p)=p;(p + \). (4.4) 

Theorem |4.2| was conjectured by Macdonald in Ref. ||M2|| (without imposing the 
Self-duality condition ( |4.3|) ). In the present self-dual set-up it is of course not nec- 
essary to distinguish between the polynomials p\(x) and the dual polynomials p\(x) 
as both polynomials coincide when Condition (|4.3|) holds. However, we have chosen 
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to keep this distinction in our notation because it is expected that with the present 
formulation all results remain valid also when the Self-duality condition ( }4.3| ) is not 
satisfied (cf. Remark [7^ of Section 0). 

Before going to the proof of Theorem |4.2| , which is relegated to Section |6], let us 



first discuss some important consequences of these duality relations. The main point 
is that the Difference equations (|3.2|) together with the Duality relations (|4.4|) imply 



a system of recurrence relations for the Koornwinder-Macdonald polynomials. To 
see this one first substitutes x = a(p + A) in the difference equations for the dual 
polynomial p*^{x), p G A. After dividing by p*^(ap) and invoking of Definition ( }4.2| ) 
one arrives at an equation of the form 

E;.(p* + fi)p;( P + \) = J2 Uj°,r-\j\(j> + x)v e j,j*(p + \)Pi(j>+\ + e e j), 

JC{l,...,n},0<|J|<r (4.5) 



ej=±l,j€J 



where 



V eJ , K {x) = YlwiEjXj) Yl ^{SjXj + EyXy)v{EjXj + EyXy + 1) 

3<i' 

X Y\ V^EjXj + Xk)v{SjXj — Xk), 



3<j' 



keK 



Uk, p ( x ) = ( _1 ) P ^2 \\ w ( £ i x i) ] [ v( £ i x i + e v x v )v{-eixi - e v x v - 1) 



L<zK,\L\=p leL IJ'eL 
ej=±l, IGL 1<V 



X 



[ | v(e t xi + x k )v(eixi- x k ), 



with 



leL 

keK\L 



viz 



shf(g* + z) 
sh(fz) 



WiZ) 



shfjg* + z) chfjgl + z) shf (g* + \ + z) chf (g* + \ + z) 
sh(fz) ch(fz) shf(± + z) chf(§ + z) 



(and E*(y) is the dual of E r (y) in Eq. (|3.2|) , i.e., with pj replaced by p*). One 
may restrict the summation in Eq. (|4.5| ) to those index sets J C {1, . . . ,n} and 
configurations of signs Sj, j G J for which A + e e j G A (|2.2|) because of the following 
lemma. 

Lemma 4.3. Let A G A = {A G Z n | A : > A 2 > ■ ■ • > A n > }. For generic 
parameters one has 

V e j,Ap + X) = iff A + e.jG'A. 
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Proof. Suppose Kj, j c (p + A) = 0. Then we are in one of the following situations. 

1. w(Ej(pj + Xj)) = for some j G J. 

For generic parameters this only happens when Xj = 0, Ej = — 1 (and j = n). 
So, A + e e j A because Xj + Ej < 0. 

2. v(ej(pj + Xj) + £j>(pj> + Xj')) = for some j,j' G J with j < j'. 

For generic parameters this only happens when —Ej = Ej> = 1, Xj = Xj> (and 
j' = j + 1). So, A + e eJ G" A because Xj + Ej < Xj/ + Ej' with j < j' . 

3. v(£j(pj + Xj) + £j'(pj' + Xj') + 1) = for some j,f G J with j < j f . 

For generic parameters this only happens when — Ej = Ej' = 1, Xj = Xj' + 1 
(and j' = j + 1). So, A + e £ j G" A because Xj + Ej < Xj> + Ej> with j < j'. 

4. v(£j(pj + Xj) + £k{pk + Afc)) = for some j G J and k G" J. 

For generic parameters this only happens when Ej = —Ek — 1, Xj — Xk and 
j — k + 1, or Ej = —£ k = —1, Xj = Xk and j — k — 1. So, A + e e j G" A because 
Aj + Ej < Xk with j E J, k <=£ J and j < A;, or Aj + Ej > A& with j G J, G" J 
and j > fc. 

Conversely, if A' = X + e £ j G" A then either A^ < or there exist aj G {1,... , n — 1} 
such that A^ < A^ +1 . The first case, i.e. A^ < 0, can occur only if n G J and 
A n = 0, E n — — 1. The vanishing of V e j,jc(p + A) then follows from the vanishing of 
yj(£ n (p n + X n )). In the second case, i.e. A^- < A'- +1 , it is convenient to distinguish the 
following three situations. 

1. jj + ieJ. 

One has A^- = Aj + Ej < Xj + i + £j + \ = Xj +1 only if < A^ — Xj + i < 1 and 
— Ej = £j + i = 1. The vanishing of V e j,j c {p + A) then follows because either 

v(Ej(pj + Xj) + Ej'{pj> + Xj')) = or V(£j(pj + Xj) + £j'{pj' + Xj') + 1) = 

(depending on whether Xj = Xj+i or A^ = Xj+i + 1). 

2. j G J, j + 1 J. 

One has A^- = Aj + Ej < Xj + i = X'j +1 only if Xj = Xj+i and Ej = —1. The 

vanishing of V e j t j c (p + A) then follows because v(Ej(pj + Xj) + (/ty + Aj')) = 

3. j g J, j + 1 G J. 

One has Aj = Aj < Aj + i + £j + i = X'j +1 only if Aj = Aj + i and Ej + i = 1. The 
vanishing of V £ j t j c (p + A) then follows because v(—(pj + Xj) + £j'(pj> + Xj')) = 0. 

(The a priori fourth situation j, j + 1 G" J does not occur because in that case 
Aj = Aj > Aj + i = Aj +1 .) □ 



In order to restrict the sum in Eq. (4J3) to the index sets and signs with A + e £ j G A 
it is of course sufficient to know that Kj, j c (p + A) = if A + e £ j ^ A. It is clear from 
the proof of Lemma |4.3| that this is actually true for all values of the parameters (the 
genericity of the parameters was needed only when proving the converse statement). 
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After restricting the sum, we may apply the duality theorem (Theorem |4.2| ) to 
obtain the equation 

E;( P * + fi)p x (p* + fi) = Yl Ujc,r~\j\(p + \)v £ji jc(p + \)p x+eEj (p* + p). 

JC{l,...,n},0<|J|<r (4.6) 
Sj=±l, jeJ; A+e eJ eA 

Equation ( |4.6| ) describes an equality between trigonometric polynomials that is sat- 
isfied for all p* + p, p G A. But then the equality must hold identically for all values 
of the argument and we arrive at the following theorem. 

Theorem 4.4 (recurrence relations). The renormalized Koornwinder-Macdonald poly- 
nomials p\(x) (f4.2| ), A G A ( |2.2| ) ; satisfy a system of recurrence relations of the form 

E*(x)p x (x)= Y Ujc >r _ lJl (p + \)V £ , lJ c(p + \)p x+eeJ (x), 

JC{l,...,n},0<|J|<r (4.7) 
£j=±l, j'SJ; A+e £j eA 

r = 1, . . . , n. 

For r = 1 the Recurrence relation fl4.7| ) becomes 
E*(x) p\(x) = 

Y w(p j +Xj)Y[ ; ij(p j +Xj+p k +X k )v(pj+X j -p k -X k ) \p x+e .(x) -p\{x)} + 

l<j<n k +3 
A+e,gA 

Y w(-pj — \j)Y\_v(—Pj-^3+Pk + *k)v(—Pj-^3—Pk — h) [Px-e^x) ~ P\(x)] 

l<j<n Mi 
A-e^eA 

(E*(y) = E*(y)). In this special case the recurrence formula was conjectured by 
Macdonald |}M2| ; after specialization to one variable (n = 1) this formula reduces to 
the three-term recurrence relation for the Askey- Wilson polynomials 



In the one- variable case both duality relations and self-duality can be easily checked 
directly through the explicit representation of pi(x) in terms of the basic hypergeo- 
metric 403-series (cf. ||AWj , Eq. (5.8)] and Relation fl2.6Q ): 

(q-i q 2 9o+ l g9o-x qgo+x \ 
_ g90+91 ' ; ?90 + S2 '+i/2 ) ^ 90+93+ i/2 ; y • ( 4 - 8 ) 

It is clear from Representation ( f4.8[ ) that the duality relations pi(g + m) = p* m (gQ + l) 
hold actually without restriction on the parameters (recall that for n = 1 one has 
p = <7q and p* = g ) and that pi(x) is self-dual iff g$ = g , or equivalently iff 
9o ~ 9i - 92 ~ 93 = (notice that g + 9r = 9o + 9* for r = 1, 2, 3). 

Remark: Recurrence relations of the type as in Theorem [4.4| are sometimes called 



generalized Pieri formulas (after similar formulas for the Schur functions) [|M3|1 . Sim- 
ilarly, the duality relations of Theorem |4.2| are also referred to as symmetry relations. 
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5. Normalization 



In this section we exploit the recurrence relations of Theorem to evaluate the 
normalization constants p\(/3p*) and (p\,P\)a (giving rise to polynomials satisfying 
the Duality relations ( |4.4j ) and orthonormal polynomials, respectively). First some 
notation is needed. Let 

A+(.x) = ] J d v (xj + Xj>) d v (xj - Xf) ] [ d w (xj), (5.1) 

with 

d ( z ) = g-9**/ 2 ( qZ > q )°° (a = e - alS ) 

(n z -n z n^/ 2+z ^ -n ( - 1 /' 2+z '>- n) 
wK ' ^ (g(ffo+ z ), -q(9t+z)^ q {g* 2 +l/2+z) ^ _ q (g*+l/2+z) . 

and let 

A + (x) = | | d v (xj + Xf) d v (xj - Xf) JJ d w (xj), (5.2) 

l<j<j'<n 

with 

w w (niz+l) fl (3/2+ Z ) _ fl (3/2+ 2 ) \ 

_ q (gZ+-+9* 3 )(z+l)/2 ' ' ' 1 ' / x 



( q (-9a+z+~L) ; g(-9 2 *+3/2+^) ) _ g (-g|+3/2+z) . 

(The function A + (x) ( |5.2| ) is obtained from A + (x) ( |5.1| ) by substituting (7* — > —(7*, 
S^-tf* (r = 0,i,2,3), andz^z + l.) 

From the relations ^(2; + 1) = {7(2;) ^(2;), + 1) = w(z)d w (z) and rf„(2 + 1) = 
v{— z — l)d v (z), d w (z + 1) = w(— z — l)d w (z) it follows that 

A+(x + e { i v .. jr} ) 



A + (x) 

and 

A + (x + e { i v .. jr} ) 



V{l,...,r},{r+l,...,n}(x) (5.3) 



V{i,...,r},{r+i,...,n}(-a; - e { i v . . jr} ), (5.4) 



A+(x) 

where Vj^x) is the same in Eq. ( |4.5| ) (with all signs Sj, j G J being positive). 
Theorem 5.1 (normalization 1: duality). One has 

f a *x A + (p + A) 

Pa(/?P ) = ~—rr — ■ (5-5) 
A +(p) 
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Proof. Expanding the l.h.s. and the r.h.s. of Recurrence relation ( |4.7|) in mono- 
mial symmetric functions and comparing the coefficients of the leading monomial 
m\ +e{l r} (x) at both sides of the equation leads to a relation between p\+ e{1 r} (flp*) 
and PxiPp*): 

VPxW) = V { i,... ,r},{r+l,...,n} 

(p + A)/p A+e{1 ,..., r} (/5/). (5.6) 
One can rewrite this relation with the aid of Eq. (5.3) in the form 

A+(p + A) A + (p + A + e {w} )' 

Because any vector in the cone A (|2.2|) can be written as a nonnegative integral 
combination of the vectors en__ j7 .y, r — 1, . . . , n, it follows from (iterated application 
of) Eq. (|5.7|) that the quotient p\((3p*)/ A + (p + A) does not depend on the choice of 
A G A. Evaluation of this quotient in A = (so p\ = 1) then entails Eq. (|5.5|). □ 



Theorem |5.1| manifestly demonstrates that p\((3p*) is indeed nonzero for generic pa- 
rameters (Lemma ^4.1|) . One observes in particular that p\((3p*) is positive for non- 
negative values of the parameters g,g r (r=0,l,2,3). (For positive parameters this is 
immediate from the expression in Eq. ( |5.5|) , whereas for g or g r equal to zero one 
needs to consider the limit g — ► + or g r — > 0+, respectively). 

Let us next turn to the computation of the (ortho)normalization constants (p\, P\) a 
(for nonnegative values of the parameters g,g , . . . , g 3 ). In the case A = (so p\ = 1), 
the corresponding (generalized Selberg-type) integral was evaluated by Gustafson in 
Ref. |GJ. In our notation, Gustafson's result may be rephrased as (cf. Remark ii., 
below) 

(l,l> A = 2 w n!A + (p)A + (p). (5.8) 

The following theorem generalizes this (constant term) formula to the case of Koorn- 
winder-Macdonald polynomials of arbitrary degree. 

Theorem 5.2 (normalization 2: orthonormality) . Let the parameters g , g r , r = 0, . . . 

be nonnegative. Then one has 

(Px, Px) a = Tn\ A + (p + A) A + (p + A). (5.9) 

Proof. Set A(x) = A + (x) A + (— x) with g*,g* replaced by g, g r (so A(x) equals A((3x), 
cf. Eq. fl2.1|) ) and let (-,-)a be the inner product determined by (cf. Eq. 



(m A ,m A /) A = — / ■•■/ m x {ix)m X '(ix) A(ix)dxi- ■ -dx n 

271 I (5.10) 



(where fh\{x) = m\(/3x)/m\(0), cf. Definition ( [4. 2D ). Clearly, the polynomials p\ 
\) are orthogonal with respect to (•, -) A (cf. Eq. ( p.5|) ). Furthermore, one has 
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(using Definition ( |4.2|) and Theorem |5.1| ) 

N x = {Px,Px)a = | f l ^, 2 (Px,Px)a = ( r A+(kP \, J (Pa> Pa) a- 

\px(Pp*)\ 2 \A + (p + \)J ( 511 ) 

By evaluating both sides of the identity 

(E*px,Px+e {1 ,..., t} )a = (P^KPX+e {1 ,..., r} ) A ( 5 ' 12 ) 

using Recurrence relation (|4.7|) and the orthogonality of the polynomials px, A G A, 
we arrive at the following relation between Nx (= (Px,Px) A ) and Nx+ e{1 r} - 

V{l,...,r},{r+l,...,n}(P + ^)Nx+e {1 r) = V{1,... ,r},{r+l,... ,n}(~ P ~ ^ ~ Hh- ,r}) N\- 

(5.13) 

With the aid of Eqs. (5.3) and (5.4) one rewrites this relation in the form 

A ± (p ± a) = A+( P + A + e , w ,) (514) 

A + (p + A) A+ *" "A + (p + A + e {1 ,... ,,)' 
which after invoking of Eq. ( |5.11| ) goes over in 

(PX, Px) A (PA+e { i,„. , r} ) PX+e (1 r} ) A 



A+(p + A) A + (p + A) A + (p + A + e { i,... , r} ) A + (p + A + e {1 ,... , r} ) ( 5 _ 15 ) 

It follows from Eq. (|5.15|) that the quotient (px,Px)a/(^+(p + A) A + (p + A)) does 
not depend on A G A fl2.2|) (cf. the proof of Theorem |5.1|) . Combining this with 
Gustafson's result (|J) for A = yields Eq. flSj). □ 



Remarks: i. Both the evaluations for px(f3p*) (Theorem |5.1|) and {px,Px)a (The- 
orem |5.2| ) were conjectured by Macdonald in Ref. ||M2|| (without imposing the Self- 



duality condition ( |4.3|) ). The formula in Theorem |5J] is sometimes referred to as the 
evaluation or specialization formula/conjecture, and the formula in Theorem |5.2| is 



also known as the inner product identity/conjecture. In the special case that A = 0, 
the latter formula is also called the constant term formula (as it amounts to explicit 
computation of the constant term in the Fourier decomposition of A(ix)). 

ii. Equation ( |5.8| ) boils down to an explicit evaluation of the generalized Selberg- 
type integral 

— ) /■■•/ A{ix)dx 1 ---dx n = 2 n n!A + (p)A + (p). (5.16) 

V27r/ J-Tr/a J-n/a 

By cancelling common factors in the numerator and denominator, it is not difficult 
to rewrite the r.h.s. of Eq. ( |5.16| ) as 

' {T,q)oc{V]q)oo[[ < r< r><3( a rar>V ,9)c 
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where t = q 9 , a = q 9 \ Oi = -q 9 \ a 2 = q 92+1 / 2 , and a 3 = ( cf Eq Q). 

This is the expression for the evaluation constant found by Gustafson Q. 

6. Proof of the duality theorem 



In this section we will prove the Duality relations ( [4.4] ) (Theorem |4.2j ) by per- 
forming induction on p. It is immediate from the definition of the (renormalized) 
Koornwinder-Macdonald polynomials that the duality relations hold for all A G 
A (|2]2|) if p = 0, as in that case both sides of Eq. ( |4.4|) are identical to one. Now, let 
w G Abe nonzero and let us assume as induction hypothesis that Eq ( [4.4|) is valid for 
all A G A and all p, G A with p < uj. We shall prove that this implies that Eq (|4.4j) 
also holds for p — uj (and all A G A). 

Since uj G A is nonzero, there must exist ansG{l,...,n — 1} such that uj s > uj s +i 

(> 0). Hence, v = uj — e{i r .. !S } G A ( [2.2|) (recall e{i v .. iS } = e± + h e s ). The most 

important step in the duality proof consists of demonstrating (of course without 
relying on Theorem |4.2| or Theorem |4.4| ) that p u (x) satisfies the Recurrence relation 
(|4.7| ) for r = s. The proof of this statement mimics the derivation presented in 
Section [| of the recurrence relations (Theorem |4.4| ) starting from the duality relations 
(Theorem [4.2| ). At some points, however, it is necessary to adapt the arguments given 



there since the starting point now is the above induction hypothesis rather than the 
duality theorem. 

As before (cf. Section f|), substituting x = a(p + u) in the s-th difference equa- 
tion for p\{x) entails (after renormalizing and restricting the sum with the aid of 
Lemma |4.^j| to those index sets and sign configurations such that v + e e j G A) 

E* s (p* + \) Pl(p + v) = Yl Ujc )a _\j ] (p + v)V eJt jc(p + v)p* x (p + v + e 6 j). 

0<\J\<s,ej=±l,jEJ (6.1) 
v+e Ej £A 

Let us assume for the moment that A < uj. Then we may use the induction hypothesis 
to rewrite Eq. (p.l|) as 



E* s {p* + A) p v {p* + A) = Ujc, a _ [J \(p + v)V B j t jc(p + v)p v+eeJ {p* + \). 

Q<\J\<s,Ej=±i,jeJ (6.2) 

i/+e eJ eA 

To obtain Eq. ( p.'2\) from Eq. (|6.1|), we have used for the l.h.s. that p* x (p+v) = p u (p* + 
A) because z/(= tu — eix^.. ,,}) < uj, and for the r.h.s. that p*x(p+v+e £ j) = p u+eeJ (p* + \) 
because either 1. A < uj, or 2. v + e £ j < uj (this happens when e £ j 7^ e{i v .. jS }), or 3. 
A = v + e £ j = uj (this happens when A = uj and e £ j = cn.^i). In the third case one 
has that p*(p + cj) = p u (p* + uj) trivially, because of the Self-duality assumption (fL3]) 
(which implies that p*(x) = pu(x) and p* = p). It is precisely at this point (and 
only at this point) that we actually use the Self-duality condition (cf. Remark |7^ of 
Section ^|). 
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So far we have shown that Eq. (|6.2| ) holds for all p* + A with A < uo. To see 
that equality actually holds for all values of the argument (and thus proving the s-th 
recurrence relation for p v {x)), we need two lemmas. 

Lemma 6.1. One has 

E*(x)p u (x) = 2J c fiPv( x ) mth c M G C. (6.3) 

Proof. It is clear from the definitions that E*(x) p v {x) is an even and permutation 
invariant trigonometric polynomial. Hence, it can be expanded as a finite linear 
combination of monomials rh\(x), A G A (where rh\(x) = m\(j3x)/m\(0), cf. Defi- 
nition ( |4 . 2| ) ) . It follows from the asymptotics at infinity that in this expansion only 



monomials rh\(x) with A < uj occur with a nonzero coefficient. To see this we set 
x = Ry with y x > y 2 > • • ■ > y n > and notice that 

m x {Ry) ~ e a(3R{x ' y) for R -> 00 (6.4) 
(where (A, y) = J2i<j< n ^jVj an< ^ ~ denotes proportionality). One furthermore has 
E*(Ry) ~ e aPR{e ^-^ y \ p v (Ry) ~ e al3R{u ' y) for R -> 00. 

(6.5) 

In obtaining the asymptotics for p u (Ry) we have used the fact that 

A' < A iff (A', y) < (A, y) G M n with yj. > • • • > y n > 0. (6.6) 

By comparing the asymptotics of E*(Ry) p u (Ry) ~ exp(aj3R(u, y)) (recall uj = v + 
e{i,..., s }) with that of fh\{Ry) (6.4), one infers — using Eq. ( |6.6|) — that 

E*(x) p v (x) = 2J c^fh^x) with c M G C. (6.7) 

The lemma is now immediate from Eq. (|6.7Q and the fact that the transition ma- 
trix relating the bases {pa^Iasa and {m x (x)} X £A is upper triangular (p\(x) = 

Ex >eA, A'<A Ca . a ' m ^( X ) m \i x ) — SA'eA,A'<A C A™A' V\'{ x ))- □ 

Lemma 6.2. For generic parameters g t go, ... ,#3 ared scale factors a, (3, the matrix 
[Pm(p* + ^)]m,Aj w ^ p, A G A ane? p, A < uj, is invertible. 

Proof. From the asymptotics fn^{p* + A) ~ exp(a/3(p, p* + A)) = g~^^*+ A ) for 
<?*(= (?) — > 00 and the fact that det[g~^ ,A ^A, M <a; 7^ for generic g (see Remark 
at the end of this section), it follows that the determinant det[m Ai (p* + \)]^\< w is 
not identically zero. The same is then true for det[p M (p* + A)] Mi a<w, as the transition 
matrix relating the bases {p\(x)}\eA an d {^a^IasA (and therefore also the matrix 
relating det[p M (p* + A)]^^ and det[m M (p* + A)] Mi a< w ) is (upper) triangular (cf. the 
proof of Lemma |6.1|). Thus, in view of the analyticity in the parameters and in the 
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scale factors a, (3, it is clear that det [p M (p* + A)]^^ 7^ generically, which proves 
the Lemma. □ 



With these two lemmas we are finally in the position to show that Eq. ( |6.2| ) not only 
holds for p* + X, X < uj but in fact for all values of the argument. Substituting x = 
p* + A in Expansion ( |6.3p and comparing with Eq. ( p.2[ ), tells us — using Lemma |6.2| — 
that the expansion coefficient c M equals Uje, a -\j\(p + u ) VeJ,j c (p + v) if p = v + e £ j 
(for some index set J with < | J\ < s and configuration of signs Ej, j G J), and that 
c M = otherwise. It is clear that with such expansion coefficients c M , Eq. (|6.3|) goes 
over in the s-th recurrence relation for p v {x): 

E*(x)p u (x) = Ujc tS _\j\(p + v)V eJi jc(p + u)p„ +eeJ (x). (6.8) 

o<|J|<s, ej=±i, jeJ 

v+e eJ <=A 

It follows in particular from the Recurrence relation ( |6.8D that Eq. (|6.2| ) is valid 
for all X £ A (and not just for A < uj). Furthermore, if we subtract Eq. (|B"72|) 
from Eq. ( |6.1| ), then all terms in the difference — except the leading term in the r.h.s. 
corresponding to the index set J = {1, . . . , s} with all signs Ej (j G J) positive— 
manifestly cancel each other in view of the induction hypothesis (recall p\{p + v) = 
p u (p* + A) because v < v + £{\,...,s) = uj, and p\(p + v + e eJ ) = p u+eeJ (p* + A) for 
e eJ 7^ e {i,... ,s} because v + e e j < v + , s } = uj if e £ j 7^ £{1,... ,«})■ Hence, we end up 
with the equation (recall C/^o = 1) 

= V {1; ... , s}j {s+li ... >n} (p + 1/) (p*(p + w) - p w (p* + A)) . (6.9) 

For generic parameters , s }, { s +i,... ,n}(p + z/ ) is nonzero (Lemma |4.3| ), so we find 
that p^(p + co>) = p w (p* + A) for all A G A (first for generic parameters and then for 
all parameters in view of the analyticity in the parameters). 



This completes the induction step and thus the proof of Theorem P . 

Remark: In the proof of Lemma |0] we have used that det[g~^'^]A iM <u; 7^ for 
generic q. This follows from a more general result due to Koornwinder | [K1| | (see also 
Ref. ||M3| , Ch. 6]) stating that for distinct vectors v\, . . . ,vn G lR n (and generic q) 

detiq^^l^N = (-lyqto'*'^ ^ 0. (6.10) 

<t£5jv 

To see that the terms q^i^^uv in Eq. (|6.10| ) cannot all cancel each other one uses the 
estimate 

^2 (Vj^)) < J2 (vj,v 3 ) (6.11) 
l<j<N i<i<v 

(immediate from the Cauchy-Schwarz inequality (twice)), with equality holding only 
when a = id. It follows from this estimate, first for q — > 00 and hence for generic 
(say positive) q by analyticity, that the determinant in Eq. ( |6.10| ) is indeed nonzero. 
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7. Concluding remarks 

7.1. Structure of the difference equations. It has already been pointed out in 
Remark i. of Section |3| that the structure of the Difference equations ( |3.2| ) is that 
of a system of eigenvalue equations for a family of commuting difference operators 
D r ( p.5| ) (r = 1, . . . , n). The heighest-order part of D r is of the form (recall Uk,o = 1) 

y i.. <■{■>■) t.,. s . (7.i) 

JC{1,... ,n}, \J\=r 
ej=±l, jeJ 

The coefficients V £ j jc{x) (|3.2| ) are related to the weight function A(x) ( |2.1|) via (cf. 
Eq. (5.3)) 

V{l,..., r },{r+l,...,n}( X ) = A + ( X + P e {h- ,r}) I A + ( X ) > ( 7 - 2 ) 

where A + (x) = A + (/3' 1 x) with replaced by g,g r (so A(a;) = A + (x)A + (—xj). 

The lower-order terms of D r (|3.5|) consist of a similar sum over the index sets J C 
{1, . . . ,n} with cardinality smaller than r and with each term V £ j j jc(x)T £ j )j 3 being 
multiplied by the function Ujc^ r _\j\(x). 

Now, for A = the r.h.s. of the Difference equation (|3.2|) vanishes because 
E r (p) = |D1|| . Equivalently, one could say that the difference operators D\, . . . , D n 
(|3.5| ) annihilate constant functions. (For D = D\ this is clear from Eq. (3.8).) This 
property of the difference operators/equations gives rise to the following set of rela- 
tions between functions V £ jj and Uk, p - 

^2 Ujc^ r _\j\ V £j ,jc = 0, r = l,...,n. (7.3) 

JC{1,... ,n},0<|J|<r 

ej=±i, jeJ 

These functional relations (together with the condition that Uk, p = 1 for p = 0) 
actually determine Uk, p completely in terms of V £ jj. For instance, for r = 1 Eq. (|7T3|) 
entails 

Ur,1 = — ^^{k}, K\{k} + V-{k},K\{k}) 
k£K 

= ~^2( W ( X k) n V(x k + X k r) V(x k - X k r) + 
keK k'eK,k'^k 

w{-x k ) Yl v(-x k +x k >)v(-x k -x k/ )j. 

k'GK, k'^k 

(Eq. ( |7.3|) immediately yields Uk,\ for K = {1, . . . the case of general K then 
follows by renumbering.) 
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More generally, one obtains Uk, p for general p from Eq. (|7.3|) with r = p by per- 
forming induction on p. This way one finds 



U 



K,p 



y] ( _1 ) m VeLi,K\Li K(L 2 \Li), AT\L 2 ' ' ' K(£ m \L m _i), l^L* 

,£J=± 

It turns out that this expression for Uk, p can be rewritten in the more compact form 
that was used in Eq. 



HlCL 1 C---CL m CK, l<m<p 

\L m \=p, £;=±1, l€Lm 



D2| , P3|| . The equality of both expressions for Uk, p hinges 
on a system of functional equations for v(z) satisfied by v(z) = sh | ((3g + z)/sh (| z). 



7.2. Dropping the self-duality condition. In Section |6| we provided a proof of the 
duality relations for the Koornwinder-Macdonald polynomials (Theorem [4.2|) with 
parameters subject to the Self-duality condition QL3Q. This condition effectively 
reduces the number of independent parameters from five to four (not counting the 
scale factors a and j3). It is expected (and conjectured by Macdonald ||M2|| ), however, 
that the Duality relations ( |4.4[ ) are true for the full five-parameter Koornwinder- 
Macdonald family (thus generalizing the state of affairs for n = 1). Should one 
succeed in proving the Duality theorem [4.2| without restrictions on the parameters, 
then automatically all other results of Sections [| and H carry over to this slightly 
more general situation (with the proofs given applying verbatim). 

A careful examination of the proof of Theorem [4.2| given in Section ^ reveals that the 
only step requiring invoking of the Self-duality condition (|4.3|) has been the derivation 
of Eq. ( |6.2|) from Eq. Q3.1|) (with the aid of the induction hypothesis). At that point 
we needed that p^(p + uj) = ftui(p* + oj), which is trivial for self-dual polynomials 
(because then p*(x) = p~u(x) and p* = p), but which requires a proof when the 
self-duality condition is dropped. If one would be able to prove the relation p* (p + 
uj) = Puj(p* + uj) (uj G A) for arbitrary parameters, then Theorem |4.2j (and thus 
all other theorems in Sections [| and ||) would follow immediately for the complete 
five-parameter Koornwinder-Macdonald family. 

There might be an alternative approach. If one chooses to define the renormalized 
Koornwinder-Macdonald polynomials p\(x) as p\(f3x) divided by the constant in the 
r.h.s. of Eq. ( |5.5| ), then the derivation of Eq. ( |6.8| ) can be established starting from 
Eq. ( |6TT| ) by assuming A < uj and applying the induction hypothesis to arrive at 
Eq. ( |6.2| ) (because A < uj it is now not necessary to invoke the self-duality condition). 
Next we bring all terms in the r.h.s. of Eq. ( |6.2| ) with | J\ — s to the l.h.s. and employ 
a version of Lemma |6.1| and Lemma with strict inequalities (p, A < uj) to arrive at 
Eq. (|6.8|). (By bringing the terms with | J\ = s to the other side one ensures that the 



resulting function in the l.h.s. can be expanded in polynomials p M (x) with p < uj.) 
This proves the induction step. However, to check now that the duality relations hold 
for p = amounts to proving Theorem |5.1| . Hence, the upshot is that all results of 
the paper can be extended to arbitrary parameters once Eq. ( |5.5|) (Theorem |5.1| ) has 
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been verified. Or, in other words, the Duality relations (|4.4|) and the Evaluation (or 



specialization) formula (5.5) follow from each other 



7.3. AfRne Hecke algebras. In Ref. ||C1|| , Cherednik proved Macdonald's orthonor- 
malization conjectures for the Macdonald polynomials related to reduced root systems 
by means of certain Hecke-algebraic techniques (leading to shift operators). More pre- 
cisely, Cherednik considered the case of admissible pairs of the form (R, R v ) with R 
reduced. Recently, these algebraic methods also resulted in a proof of the correspond- 
ing duality relations and renormalization formulas (i.e., the evaluation/specialization 
formulas, cf. Remark i. of Section |5|.) ||C2|| . For the type A root system this approach 
reproduces (albeit in a completely different manner) the results of Koornwinder and 
Macdonald ||K1|, |M3|| . For the remaining classical reduced root systems (i.e., the 
types B, C, and D), the (R, i? v )-type Macdonald polynomials amount to (self-dual) 
Koornwinder-Macdonald polynomials with special parameters g r (type B: g = g 2 , 
91,92, = 0; type C: g = g x , g 2 , g 3 = 0; type D: g Q , g u g 2 , g 3 = 0). (For the root 
systems B n and D n one obtains only half the polynomials via the above specializa- 
tion of the parameters; the other half of the polynomials is obtained by specializing 
to parameters that are not self-dual, cf. Ref. |pi| , Section 5].) It hence follows that 



for the above parameters Refs. [CT| , |C2|| provide an alternative approach towards the 
proof of the Theorems [4.2| , |5.1| , and |5.2| . 

Meanwhile, both Noumi |[N2j| and Macdonald |[M4j| independently announced that 
Cherednik's Hecke-algebraic techniques may be extended to the full five-parameter 
Koornwinder-Macdonald family. Hopefully this will eventually lead to an alternative 



proof of the Theorems [4.2| , |5.1| , and |5.2| valid for all parameters. 

One aspect of the present paper seems hard to achieve via Hecke algebras, though. 
It does not seem very likely that such algebraic methods will independently reproduce 
our explicit formulas for the difference equations and the recurrence relations for 
the Koornwinder-Macdonald polynomials. At present Cherednik's techniques permit 
concluding the existence of such difference equations and recurrence relations (of 
course for special parameters) without providing much clue as regards to their explicit 
form except in the simplest cases, e.g. when r = 1 (cf. in this connection also the 
statements in second paragraph of the proof of Theorem 5.1 of Cherednik's paper 
Ref. £3). 

7.4. Quantum groups. Another algebraic framework in which the Koornwinder- 
Macdonald polynomials appear is the representation theory of (compact) quantum 
groups. Specifically, it has been demonstrated recently by Noumi and Sugitani that 
for certain special values of the parameters the Koornwinder-Macdonald polynomi- 
als may be interpreted as zonal spherical functions on compact quantum symmetric 
spaces of classical type ||NSj1 (see also Ref. [[Nil for a detailed treatment of the type 
A case). 
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In addition, it turned out that for the root system A n _i Macdonald's polynomials 
may also be realized as vector valued characters of U q (sl n ) [|EK1 ]. This observation 
has led to yet another (representation-theoretic) proof of Koornwinder's duality and 
recurrence relations for the A-type Macdonald polynomials ||EK2| . 



7.5. Integrable systems. It is possible to view the commuting difference operators 
Dx, . . . , D n fl3.5| ) as a complete set of quantum integrals for an integrable quantum 
mechanical n-particle model | |D2| , |D3| |. Similar integrable systems associated with 
and diagonalized by the Macdonald polynomials related to classical root systems are 
obtained via limit transitions (type A) or specialization of the parameters (type B, C, 
D, and BC) ||L)1| - For the type A root systems the commuting difference operators 
(quantum integrals) of the model were already found independently by Ruijsenaars 



Q and Macdonald |Ml[ 



In this special case, Ruijsenaars also studied in great detail 

It is interesting 



the properties of the corresponding classical mechanical systems |R2 



to note that also at the classical level duality relations, which were actually known 
even before their quantum counterparts were discovered, play a crucial role in solving 
the system. 
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